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I. INTRODUCTION 
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We investigate gravitational lensing in the Palatini approach to the f(R) extended theories of 
gravity. Starting from an exact solution of the f(R) field equations, which corresponds to the 
Schwarzschild-de Sitter metric and, on the basis of recent studies on this metric, we focus on some 
lensing observables, in order to evaluate the effects of the non linearity of the gravity Lagrangian. 
We give estimates for some astrophysical events, and show that these effects are tiny for galactic 
lenses, but become interesting for extragalactic ones. 
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It is well known and accepted nowadays that the Universe is undergoing a phase of accelerated expansion [l[. 
General Relativity (GR), in its standard formulation, cannot explain this fact without requiring the existence of the 
O ,. so called dark energy, or introducing a cosmological constant which, in turn, brings about other problems, concerning 
its nature and origin Q. According to a different viewpoint, current observations might suggest the failure of GR to 
describe gravitational interactions at cosmological scale: as a natural consequence, the claim arises for a modification 
of the theory of gravity. In recent years, there has been much interest for the so called f(R) theories of gravity: in 
these theories, the gravitational Lagrangian depends on an arbitrary analytic function / of the scalar curvature R 
(see Q and references therein) . These theories are also referred to as extended theories of gravity, since they naturally 
| generalize GR on a geometric ground: in fact, when f(R) = R the action reduces to the usual Einstein-Hilbert action, 
and Einstein's theory is obtained. The field equations of these theories can be obtained in the metric formalism [J], by 
varying the action with respect to metric tensor, or in the Palatini formalism, where the action is varied with respect 
t-H , to the metric and the affine connection, which are supposed to be independent from each other [5j. The increasing 
interest for these theories is mainly due to the fact that it seems that they are able to provide cosmologically viable 
models (see @ and references therein), in which both the late time acceleration and the inflation phase are present. 
However, we must remember that even though there are problems in explaining the cosmological dynamics, GR 
t-H ' is in excellent agreement with other gravity tests, such as Solar System ones and those coming from binary pulsar 
observations 0]; hence, as a matter of fact, every modified theory of gravity should reproduce GR in a suitable weak 
field limit. This issue is crucial for f(R) theories: their weak field limits were largely analyzed, and the comparison 
with some known tests of gravity was carried out. A lively debated developed because some models have been ruled 
out H, while other seem to pass Solar System tests jsj. Among the other things, stellar models have been investigated, 
and there is not a common agreement on the results. In fact, on the one hand, the Palatini approach seems to provide 
■ self consistent solutions llCfl. on the other hand, the same problem studied in the metric approach leads to unacceptable 
results [ll[ (however see|12j: see [l3[ for general difficulties with the Palatini approach). Gravitational lensing (i.e. 
the deflection of light rays due to the curvature of space-time), besides having an historical importance, since it gave 
one of the first confirmations of Einstein's theory in the famous eclipse of 1919, is very important today, because data 
coming from modern lensing observations are crucial to explore the Universe at very different scales [3, EE EH . It is 
then natural to ask whether modifications of the gravity Lagrangian, as those introduced in the f{R) theories, affect 
the standard theory of lensing (where "standard" stands for "GR lensing") and, in particular, the interpretation of 
its phenomenology. Gravitational lensing in the framework of f(R) was studied in |17| : in particular, theories with 
f(R) oc R n were investigated in the metric approach, with emphasis on point-like lenses. In the present paper we 
want to address a similar task, but working in the Palatini approach (however, as we remark below, our results apply 
also to a subset of the solutions metric f(R) gravity). In a previous work |18j | we found an exact solution of the 
f(R) field equations for the vacuum case, which corresponds to the Schwarzschild-de Sitter metric: starting from 
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this solution and taking into account recent studies on this metric [19j . we want to investigate the impact of the non 
linearity of the gravity Lagrangian on some lensing observables. The paper is organized as follows: in Section |TT] we 
briefly review the Palatini approach to f(R) theories, which lead to the Schwarzschild-de Sitter vacuum solution; in 
Section HTT1 we focus on point like-lensing, in GR (Schwarzschild space-time) and in f(R) theories (Schwarzschild-de 
Sitter space-time). Finally discussions and conclusions are outlined in Section HVl 



II. A VACUUM EXACT SOLUTION OF THE f(R) GRAVITY FIELD EQUATIONS 

The equations of motion of f(R) gravity in the Palatini formalism can be obtained by independent variations with 
respect to the metric and the connection from the action 1 

A = A giav + A mat = J [y/gf(R) + 2 K £ mat (<A, V^)] d 4 z, (1) 

where R = R(g,T) — g a @ R a p(T), i?^„(r) is the Ricci-like tensor of any torsionless connection T independent from 
the metric g, which is assumed here to be the physical metric. The gravitational part of the Lagrangian is represented 
by any real analytic function f(R) of the scalar curvature R. The matter Lagrangian L ma t is functionally depending 
on unspecified matter fields \& together with their first derivatives, equipped with a gravitational coupling constant 



/■>■ = 87r 4 G . 

According to the Palatini formalism [5], from (fl]) we obtain the following equations of motion: 

f{R)R { ^ ) {T)- l -f{R)g^ = kT™*, (2) 
V r a [^gf'(R)gn = 0, (3) 

where T^ at — — -^= ^ mat denotes the matter source stress-energy tensor and V r means covariant derivative with 

respect to the connection T. The equations of motion ((2|) can be supplemented by the scalar- valued equation obtained 
by taking the g-trace of ©, where we set t = trT = g^TJ^f 1 : 

f{R)R-2f{R) = KT. (4) 

The algebraic equation is called the structural equation and it controls the solutions of equations 

The field equations ([2][3]) and the structural equation ([4]) in vacuum become 

[f'(R)}R ifiu) (T)-^[f(R)]g^ = 0, (5) 

<(V9 lf(R)} gn - o, (6) 

f(R)R-2f(R) - 0. (7) 
As shown in [l8[ (see, in particular, Section 3), the system of equations ([5][7]) has the spherical symmetrical solution 

clr 2 

ds 2 = -b(r)dt 2 + -TT-+ r 2 dtf 2 + r 2 sin 2 dd<p 2 , (8) 
b(r) 

where b(r) = (l — ^jf — ; M is the mass of the source of the gravitational field and k = —Ci/4, where R = ci 

is any of the solutions of the structural equation ([7]). In doing so, we have obtained a solution with constant scalar 
curvature R. In particular, if f(R) = R (i.e. our theory is GR) then R = is the solution of the structural equation, 
and ([S]) reduces to the classical Schwarzschild solution. 

Remark. We would like to point out that, for a given f(R) function, in vacuum case the solutions of the field 
equations of Palatini f(R) gravity are a subset of the solutions of the field equations of metric f(R) gravity [2(J. In 
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particular the spherically symmetric solution ([8]) is also a solution of the field equations of metric f(R) gravity with 
constant scalar curvature R: actually, in metric f(R) gravity there are also spherically symmetric solutions where 
the scalar curvature is not constant (see e.g. (2l| ). So, strictly speaking, what follows applies also to metric f(R) 
gravity. 

In GR, the Schwarzschild-de Sitter (JSj> solution corresponds to a spherically symmetric solution of Einstein field 
equations, with a cosmological term Ag M „, A being the cosmological constant. In practice, it is A = k in our notation. 
As we said, the cosmological constant is one of the candidates for explaining the accelerated expansion of the Universe 
(in particular, a positive cosmological constant is required). The data suggest that it cannot exceed the upper value 
of A ~ 10~ 52 ?7i~ 2 [22]. As a consequence, it is reasonable to estimate that the solutions of the structural equation 
have the same order of magnitude; the same holds for the parameter k in the metric @. 

In a recent paper (l9j it has been shown that, contrary to previous claims (e.g . J23I) , the k terms in ||SJ) does 
contribute to the observed bending of light rays (see also the subsequent papers [2J], [25|] , [26| , [13] ) • The key point is 
that the Schwarzschild-de Sitter space-time does not become flat at spatial infinity and, then, the spatial geometry 
affects deflection measurements. 

In our approach, the parameter k is ultimately related to the non linearity of the action, in fact when f(R) = R — > 
k = 0. Hence, the study of the effects of the k term on light bending can be interpreted as an investigation of the 
effects of the non linearity of the gravity Lagrangian, which is what we want to address here. 

According to fl9( j . the bending angle of a light ray in the gravitational field of a massive object described by ©, 
turns out to be 

up to first order in both M and k, if both the source of light rays and the observer are far from the massive object 
(the lens). We see that, when k — 0, i.e. in GR, esds reduces to the classical bending angle egr = ^jf (e.g. [HI). The 
parameter p, to lowest order in M and k is equivalent to the impact parameter. In what follows, we will investigate 
gravitational lensing on the basis of the bending angle © . 

III. LENSING FROM POINT-LIKE SOURCES 

In this Section, we start by briefly reviewing the basic notions of gravitational lensing from a spherically symmetrical 
point-like source in GR, i.e. we consider lensing in the Schwarzschild field (for a thorough study of Schwarzschild 
lensing, see e.g. [2!|) Then, we consider the case of f(R) theories which, on the basis of what we have said before, 
corresponds to studying lensing in the Schwarzschild-de Sitter space-time. 

A. Basics of lensing from a spherically symmetrical source 

In Figure [1] a simplified lensing geometry is depicted. The observer in O sees the image of the source, located at 
S, as if it were in Si. The lens is located at the point L. In the so called thin lens approximation (see [Til . [TH [iH ]) 
the light paths are approximated by straight lines; if the lens has spherical symmetry, which is what we assume here, 
light rays propagate in a plane. We may write the following relation (see e.g. (29[), which links the angular positions 
of the images ($), the actual position of the source (j3) and the bending angle (e) 

D s t&itd = D s ta.it [3 + D ds terns, (10) 

and, since D s tan a = Dd s tan e: 

Z3 s tani9 = -D s tan/3 + -D s tana. (11) 
So, if we assume that the angles are small (weak lens approximation), we are able to write the following lens equation 



(12) 
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FIG. 1: A simplified lensing geometry: light paths are depicted as straight lines. The observer in O sees the image of the 
source, located at S, as if it were in Si. The lens is located at the point L. Distances between the observer and the source, the 
observer and the lens and the lens and the source are, respectively, defined by D s , Dd, Dd s - 



The latter equation allows to obtain the angular positions of the images once those of the sources and the bending 
angle are known. On using back the fact that aD s ~ eDd s , we may write 

= 0- ^e(0), (13) 

in which the deflection angle s appears. 

For instance, in the case of lensing in the Schwarzschild gravitational field 



dr 2 

ds 2 = -a{r)dt 2 + -7- + r 2 dd 2 + r 2 sin 2 tfdip 2 , (14) 
a{r) 

where a(r) = (l — The point-like lens is a spherical object of mass M; the bending angle turns out to be 

e = ecu = 4M//J, where p = Dd'd is the impact parameter. On substituting, we get the explicit form of the lens 
equation in this case 

The position corresponding to j3 = 0, which describes the alignment of observer, lens and source, defines the so 
called Einstein angle 



= V 4M lt? <16) 

which allows to define the characteristic length scale of lensing 

Re = # E D d , (17) 
which is called Einstein radius. On introducing the Einstein angle in the lens equation l|15p. we get 

f3tf = $ 2 --d 2 E , (18) 

or expressing angles in units of the Einstein angle #e, i.e. setting (3 = ^^s,t9 = W^, the lens equation has the 
following expression 

/M = 2 -1. (19) 
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The latter equation can be simply solved to give the angular positions of the images 



K = (20) 
and the angular separation between the two images 

Atf = ■&+ - = \Jj3 2 +4. (21) 

B. Lensing in the Schwarzschild-de Sitter field 

In order to obtain the explicit expression of the lens equation in the Schwarzschild-de Sitter space-time, we have to 
introduce the bending angle ((9} in eq. (fT3|) . So, taking into account that p = Dd$, we may write 

n „ D ds AM D ds D3$ 3 k 

^*-DjT d — + ^r L wr (22) 



On introducing the following dimensionless parameter 

2 D 2 n 2 

X ^ f dT' (23) 

and expressing angles in units of the Einstein angle $ e defined by (|16p , the lens equation becomes 

fiti = d 2 - 1 + d 4 X - (24) 

The lens equation differs from the GR one (TIT))) because of the term proportional to \ which, in turn, is 
proportional to the parameter k, related to the non linearity of the gravity Lagrangian. Indeed, we notice that, 
beyond the parameter k, x depends on the geometry of lensing, through the distances D ds , D d , D s , which enter its 
definition; \ is reasonably small in actual astrophysical events, chiefly because of the smallness of the parameter k, 
whose expected order of magnitude is the one of the cosmological constant (see Section |TT|. As a consequence, the 
term $ 4 x can be treated as a perturbation with respect to the GR lens equation (fl"9j) . 

In what follows, we are going to analyze some obscrvables in gravitational lensing, in order to evaluate whether the 
non linearity of the Lagrangian (TTJ) raises detectable effects. 

To begin with, we calculate the solutions of the lens equation (|24]) for the case = 0, which corresponds to perfect 
alignment of observer, lens and source (superior conjunction): in other words, we look for the modified Einstein angle 
in f(R) theories. 

In what follows we use "&e,gr to refer to the standard Einstein angle (|16[) . i.e. the one obtained in GR, for a 
spherically symmetrical point-like lens (described by the gravitational field ifM]) ), while we use i?e,etg to indicate 
the Einstein angle in f(R) for a spherically symmetrical point-like lens (described by the field 1(8)1). Furthermore, we 
write angles in units of $e,gr = #e ; defined in ((To]) . 

We look for a solution of the equation 

= ¥ - 1 + d\, (25) 

in the form 

$E,ETG = 1 + 6, (26) 

where S is assumed to be small, i.e. of first order in \- We then obtain 

* = (27) 

Hence we have, up to first order in \ 

$e,etg — $e,gr _x 
Ve,gr ~ 2" 



(28) 



In particular, we see from ((2~7| that the modified Einstein angle , 0e,etg is smaller (bigger) than the classical 
one "&e,gr if X is positive (negative): this means that, according to (|23[) . a positive (negative) value of k decreases 
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(increases) the Einstein angle. If we recall that k can be interpreted as a cosmological constant (see Section ITT1) . 
and we remember that it is required to be positive to mimic the accelerated expansion of the universe, we conclude 
that in this approach to f(R) we obtain an Einstein angle that is smaller than the standard GR one: this is indeed 
independent on the analytic form of the f{R), and depends on the solutions of the structural equation only. 

We continue our analysis by taking into account the lensing-induced magnification: actually, gravitational lenses 
can amplify the luminosity of the sources, so that faint sources can become visible. In the case of point-like lenses, 
the magnification turns out to be [l5| 



A 



(3df3 



(29) 



Then, from the lens equation (|24| we obtain, up to first order in x 



(30) 



We see that when \ = the GR result A = 1 1 — ^ I is recovered. Actually, in observations what is seen is the 
total luminosity, i.e. the sum of the magnification of each image. However, in practice, it is not always possible to 
see the (multiple) images of the sources, together with their magnification, but the magnification can be detected 
thanks to the relative motion of the lens and the source, that gives rise to a lensing-induced time variability in the 
magnification. The time scaling for these variations can be estimated by tE — Re/v = D^e/v, where v is a typical 
velocity of the lens (transverse to the line of sight). Consequently, a variation of the Einstein angle provokes a variation 
of this time scale. If we set tETG — Dci^etg/ v ^ according to what we have seen, then, we can estimate the variation 
as 

M = tETC-tE =5= _X_ (31) 
tE tE 2 

Let us now focus on the estimate of the lens mass, determined by lensing phenomenology: if we suppose that the 
distances and the transverse velocity are known, a measurement of tE gives an estimate of the mass M of the lens. 
As a consequence if we "incorrectly" use the Einstein angle of GR instead of the f(R) one, we get a different (biased) 
mass of the lens. If we denote as M the "true" mas of the lens, and as Mqr the one obtained by using the GR 
expression of the Einstein angle, we get 

So, according to what we have said above, the effects of f(R) is an underestimate of the lens mass, in other words 
the true lens mass is bigger than the one estimated in standard GR lensing. 

Another phenomenological issue that can be investigated is how the positions of the images (|20[) are modified in 
f(R) theories of gravitation. Since \ is reasonably small, the lens equation (|24p can be solved perturbatively: in other 
words we look for the lowest order corrections of the k term on the GR solutions. To this end, let $o be a solution of 
the lens equation corresponding to \ = 0: 

0§ o = 0§ - 1. (33) 
In particular, 0o is equal to the solutions (|20p given above. Then, we look for solutions in the form 

= o + 0i, (34) 
where 0i is supposed to be small. On substituting (|34)) in ([24"]) . and linearizing in 0i, we obtain, up to first order in \ 

01 = J- (35) 



As a consequence, we may write 



0± = ^o± - - 0± - x- (36) 
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D d [Kpc] D a [Kpc] D ds [Kpc] 


X 


20 


50 


30 


8.6 x 10" 


10 


7.6 


7.6 


10 -3 


6.0 x 10" 


18 


7.6 


50 


40 


2.5 x 10" 


10 


5 x 10 3 


10 4 


10 4 


1.5 x 10" 


-4 



TABLE I: Evaluation of the dimensionless parameter \ (which parameterizes the deviation of the f(R) theories from GR) 
as a function of the distances between the observer and the lens, the observer and the source and the lens and the source, 
respectively defined by Dd, D s , Dd s (all measured in Kpc); the value of the parameter k is 10~ 52 m~ 2 . 



If we set A$o = V fl 2 + 4 according to (|2"Tj) . we may write the angular separation between the two images in the 
form 

Atf = Ado - s(Atf )x, (37) 

where 



?(Atf ) = 



Az9§ - 4 



Atf 



1 



16A?9o 



(38) 



In our approach, the modifications Q28p. (|3"0)) , (f3"Tj) , (f3"2")> . ([57)) determined by the non linearity of the gravity 
Lagrangian consist in terms that are proportional to the parameter x, which is what one would expect in a first order 
approximation. As a consequence, the evaluation of the detectability of these modifications requires an estimate of 
the parameter x m some actual astrophysical events, which will be done in the following Section. 



IV. DISCUSSION AND CONCLUSIONS 



We have evaluated some effects of the non linearity of the gravity Lagrangian on lensing phenomenology, in the 
Palatini approach to f(R) extended theories of gravity, however our results apply also to a subset of solutions of 
metric f(R) gravity. 

These effects are controlled by the parameter Xj the latter is proportional to k which, in turn, is simply related to 
the solution of the structural equation ([JJ and it parameterizes the deviation of the f(R) theories from GR: namely, 
when k = 0, GR is recovered. The parameter x depends also on the lens geometry, through the distances Dd s , Dd, D s , 
which enter its definition (f2"3")) ; x is expected to be small, because of the smallness of the parameter k, whose order 
of magnitude is comparable to the cosmological constant. Taking into account this fact, we have evaluated the 
f(R) effects as perturbations of the standard GR results. To quantitatively evaluate these effects it is necessary to 
give numerical estimates of the dimensionless parameter x m some astrophysical events. Taking into account some 
reasonable values for lenses in our galaxy, we obtain numerically: 



To fix the ideas, we can analyze some specific events. To begin with, let us consider a lens in the galactic halo, with 
Dd = 20 Kpc, and and extra galactic source at D s — 50 Kpc: we obtain x = 8.6 x 10 -10 . Then, we can consider the 
supermassive black hole in the radio source Sgr A*, with Dd — 7.6 Kpc, and a source (1) just behind the black hole 
D s = 1 pc, (2) outside the galaxy D s = 50 Kpc we get, respectively: (1) x = 6.0 x 10~ 18 , (2) x = 2.5 x lO" 10 . 

In general, for sources outside the galaxy we may set D s ~ Dd s ; consequently if we think of distant galaxies (or 
cluster of galaxies) acting as lenses, we obtain numerically 
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For instance if we consider Dd = 5 Mpc, we obtain % = 1.5 x 10 -4 . These estimates are summarized in Table Q] 
We remember that the order of magnitude of the lcnsing phenomena is roughly determined by the Einstein angle; 
according to eq. ([28| modifications of the Einstein angle due to the non linearity of the gravity Lagrangian are 
proportional to X®e,gr- As a consequence, the smallness of \ makes the effects of f(R) theories on lcnsing hardly 
detectable for galactic lenses. The situation is different for extragalactic lenses: in fact, in this case we have \0e.gb oc 

M 1 / 2 !)^ 2 , so that distant galaxies or clusters of galaxies acting as lenses can be interesting for our purposes. We notice 
that similar results have been obtained in [24j where, working in the GR framework, the effects of the cosmological 
constant on lensing by distant clusters of galaxies have been evaluated. However, in the case of extragalactic sources, 
i.e. when the lenses are galaxies or clusters of galaxies, the point-like lens model that we have used in this paper might 
be oversimplified: a more realistic model is needed in order to investigate the possibility of detecting these effects by 
means of the present and foreseeable observational techniques. 

Nonetheless the estimates that we have obtained suggest at least the scale where the impact of the /(i?)-induccd 
corrections on lensing phenomenology becomes interesting. 

To summarize, the non linear terms in the gravity action determine a modification of lensing phenomenology; 
numerically for actual astrophysical events, these modifications are very small if we deal with lenses in our galaxy: 
these results suggest that the effect of f{R) are confined around a cosmological scale and, hence, they are not 
effective at the galactic scale. On the contrary, in our simple model, extragalactic lenses lend themselves as natural 
candidates for studying the effects of extended theories of gravity. However, more realistic models and further studies 
are needed to evaluate the detectability of these effects. 
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